Axial precession (astronomy)

For other types of astronomical precession, see Precession: Astronomy
Precessional movement.

In astronomy, axial precession is a gravity-induced, slow and continuous change in the orientation of an astronomical body's rotational axis. In particular, it refers to the gradual shift in the orientation of Earth's axis of rotation, which, like a wobbling top, traces out a cone in a cycle of approximately 26,000 years (called a Great or Platonic year in astrology). The term "precession" typically refers only to this largest secular motion; other changes in the alignment of Earth's axis — nutation and polar motion — are much smaller in magnitude.

Earth's precession was historically called precession of the equinoxes because the equinoxes moved westward along the ecliptic relative to the fixed stars, opposite to the motion of the Sun along the ecliptic. This term is still used in non-technical discussions, that is, when detailed mathematics are absent. Historically,[1] Hipparchus is credited with discovering precession of the equinoxes. The exact dates of his life are not known, but astronomical observations attributed to him by Ptolemy date from 147 BC to 127 BC.

With improvements in the ability to calculate the gravitational force between planets during the first half of the 19th century, it was recognized that the ecliptic itself moved slightly, which was named planetary precession as early as 1863, while the dominant component was named lunisolar precession.[2] Their combination was named general precession instead of precession of the equinoxes. Lunisolar precession is caused by the gravitational forces of the Moon and Sun on Earth's equatorial bulge, causing Earth's axis to move with respect to inertial space. Planetary precession (actually an advance) is due to the small angle between the gravitational force of the other planets on Earth and its orbital plane (the ecliptic), causing the plane of the ecliptic to shift slightly relative to inertial space. Lunisolar precession is about 500 times larger than planetary precession.[3] In addition to the Moon and Sun, the other planets also cause a small movement of Earth's axis in inertial space, making the contrast in the terms lunisolar versus planetary misleading, so in 2006 the International Astronomical Union recommended that the dominant component be renamed the precession of the equator and the minor component be renamed precession of the ecliptic, but their combination is still named general precession.[4]

Contents

Effects

The precession of the Earth's axis has a number of observable effects. First, the positions of the south and north celestial poles appear to move in circles against the space-fixed backdrop of stars, completing one circuit in 25,771.5 years (2000 rate). Thus, while today the star Polaris lies approximately at the north celestial pole, this will change over time, and other stars will become the "north star".[1] As the celestial poles shift, there is a corresponding gradual shift in the apparent orientation of the whole star field, as viewed from a particular position on Earth.

Secondly, the position of the Earth in its orbit around the Sun at the solstices, equinoxes, or other time defined relative to the seasons, slowly changes.[1] For example, suppose that the Earth's orbital position is marked at the summer solstice, when the Earth's axial tilt is pointing directly towards the Sun. One full orbit later, when the Sun has returned to the same apparent position relative to the background stars, the Earth's axial tilt is not now directly towards the Sun: because of the effects of precession, it is a little way "beyond" this. In other words, the solstice occurred a little earlier in the orbit. Thus, the tropical year, measuring the cycle of seasons (for example, the time from solstice to solstice, or equinox to equinox), is about 20 minutes shorter than the sidereal year, which is measured by the Sun's apparent position relative to the stars. Note that 20 minutes per year is approximately equivalent to one year per 25,771.5 years, so after one full cycle of 25,771.5 years the positions of the seasons relative to the orbit are "back where they started". (In actuality, other effects also slowly change the shape and orientation of the Earth's orbit, and these, in combination with precession, create various cycles of differing periods; see also Milankovitch cycles. The magnitude of the Earth's tilt, as opposed to merely its orientation, also changes slowly over time, but this effect is not attributed directly to precession.)

For identical reasons, the apparent position of the Sun relative to the backdrop of the stars at some seasonally fixed time, say the vernal equinox, slowly regresses a full 360° through all twelve traditional constellations of the zodiac, at the rate of about 50.3 seconds of arc per year (approximately 360 degrees divided by 25,771.5), or 1 degree every 71.6 years.

For further details, see Changing pole stars and Polar shift and equinoxes shift, below.

History

Hellenistic world

Hipparchus

Though there is still-controversial evidence that Aristarchus of Samos possessed distinct values for the sidereal and tropical years as early as ca. 280 BC,[5] the discovery of precession is usually attributed to Hipparchus (190–120 BC) of Rhodes or Nicaea, a Greek astronomer. According to Ptolemy's Almagest, Hipparchus measured the longitude of Spica and other bright stars. Comparing his measurements with data from his predecessors, Timocharis (320–260 BC) and Aristillus (~280 BC), he concluded that Spica had moved 2° relative to the autumnal equinox. He also compared the lengths of the tropical year (the time it takes the Sun to return to an equinox) and the sidereal year (the time it takes the Sun to return to a fixed star), and found a slight discrepancy. Hipparchus concluded that the equinoxes were moving ("precessing") through the zodiac, and that the rate of precession was not less than 1° in a century, in other words completing a full cycle in no more than 36000 years.

Virtually all Hipparchus' writings are lost, including his work on precession. They are mentioned by Ptolemy, who explains precession as the rotation of the celestial sphere around a motionless Earth. It is reasonable to assume that Hipparchus, like Ptolemy, thought of precession in geocentric terms as a motion of the heavens.

Ptolemy

The first astronomer known to have continued Hipparchus' work on precession is Ptolemy in the 2nd century. Ptolemy measured the longitudes of Regulus, Spica, and other bright stars with a variation of Hipparchus' lunar method that did not require eclipses. Before sunset, he measured the longitudinal arc separating the Moon from the Sun. Then, after sunset, he measured the arc from the Moon to the star. He used Hipparchus' model to calculate the Sun's longitude, and made corrections for the Moon's motion and its parallax (Evans 1998, pp. 251-255). Ptolemy compared his own observations with those made by Hipparchus, Menelaus of Alexandria, Timocharis, and Agrippa. He found that between Hipparchus' time and his own (about 265 years), the stars had moved 2°40', or 1° in 100 years (36" per year; the rate accepted today is about 50" per year or 1° in 72 years). He also confirmed that precession affected all fixed stars, not just those near the ecliptic, and his cycle had same period of 36000 years as found by Hipparchus.

Other authors

Most ancient authors did not mention precession and perhaps did not know of it. Besides Ptolemy, the list includes Proclus, who rejected precession, and Theon of Alexandria, a commentator on Ptolemy in the 4th century, who accepted Ptolemy's explanation. Theon also reports an alternate theory:

According to certain opinions ancient astrologers believe that from a certain epoch the solstitial signs have a motion of 8° in the order of the signs, after which they go back the same amount. . . . (Dreyer 1958, p. 204)

Instead of proceeding through the entire sequence of the zodiac, the equinoxes "trepidated" back and forth over an arc of 8°. The theory of trepidation is presented by Theon as an alternative to precession.

Alternative discovery theories

Babylonians

Various claims have been made that other cultures discovered precession independent of Hipparchus. At one point it was suggested that the Babylonians may have known about precession. According to Al-Battani, the Chaldean astronomers had distinguished the tropical and sidereal year (the value of precession is equivalent to the difference between the tropical and sidereal years). He stated that they had, around 330 BC, an estimation for the length of the sidereal year to be SK = 365 days 6 hours 11 min (= 365.258 days) with an error of (about) 2 min. It was claimed by P. Schnabel in 1923 that Kidinnu theorized about precession in 315 BC. Otto Neugebauer's work on this issue in the 1950s superseded Schnabel's (and earlier, Kugler's) theory of a Babylonian discoverer of precession.[6]

In recent decades, the hypothesis was revived and amplified in de Santillana and von Dechend's Hamlet's Mill (Harvard University Press, 1969). In an application of extreme Panbabylonism to archaeoastronomy, they proposed that a Babylonian mythological account of the precession gave rise via diffusion to similar myths around the world, even as far away as China, Polynesia, and North America. While their theory has not been widely accepted in academia, it anticipated the recent popular revival of interest in precessional archeoastronomy.

Ancient Egyptians

Similar claims have been made that precession was known in Ancient Egypt prior to the time of Hipparchus, but these remain controversial. Some buildings in the Karnak temple complex, for instance, were allegedly oriented towards the point on the horizon where certain stars rose or set at key times of the year. A few centuries later, when precession made the orientations obsolete, the temples would be rebuilt. However, the observation that a stellar alignment has grown wrong does not mean that the Egyptians understood that the stars moved across the sky at the rate of about one degree per 72 years. Nonetheless, they kept accurate calendars and if they recorded the date of the temple reconstructions it would be a fairly simple matter to plot the rough precession rate. The Dendera Zodiac, a star-map from the Hathor temple at Dendera from a late (Ptolemaic) age, supposedly records precession of the equinoxes (Tompkins 1971). In any case, if the ancient Egyptians knew of precession, their knowledge is not recorded in surviving astronomical texts.

Michael Rice wrote in his Egypt's Legacy, "Whether or not the ancients knew of the mechanics of the Precession before its definition by Hipparchos the Bithynian in the second century BC is uncertain, but as dedicated watchers of the night sky they could not fail to be aware of its effects." (p. 128) Rice believes that "the Precession is fundamental to an understanding of what powered the development of Egypt" (p. 10), to the extent that "in a sense Egypt as a nation-state and the king of Egypt as a living god are the products of the realisation by the Egyptians of the astronomical changes effected by the immense apparent movement of the heavenly bodies which the Precession implies." (p. 56) Following Carl Gustav Jung, Rice says that "the evidence that the most refined astronomical observation was practised in Egypt in the third millennium BC (and probably even before that date) is clear from the precision with which the Pyramids at Giza are aligned to the cardinal points, a precision which could only have been achieved by their alignment with the stars. This fact alone makes Jung's belief in the Egyptians' knowledge of the Precession a good deal less speculative than once it seemed." (p. 31) The Egyptians also, says Rice, were "to alter the orientation of a temple when the star on whose position it had originally been set moved its position as a consequence of the Precession, something which seems to have happened several times during the New Kingdom." (p. 170) see also Royal Arch and the Precession of the Equinoxes

The notion that an ancient Egyptian priestly elite tracked the precessional cycle over many thousands of years plays a central role in the theories expounded by Robert Bauval and Graham Hancock in their 1996 book Keeper of Genesis. The authors claim that the ancient Egyptians' monumental building projects functioned as a map of the heavens, and that associated rituals were an elaborate earthly acting-out of celestial events. In particular, the rituals symbolised the "turning back" of the precessional cycle to a remote ancestral time known as Zep Tepi ("first time") which, the authors calculate, dates to around 10,500 BC.

Mayans

There has been speculation that the Mesoamerican Long Count calendar is somehow calibrated against the precession, but this view is not held by professional scholars of Mayan civilization.[7]

Indian views

A twelfth century text by Bhāskar II[8] says: "sampāt revolves negatively 30000 times in a Kalpa of 4320 million years according to Suryasiddhanta, while Munjāla and others say ayana moves forward 199669 in a Kalpa, and one should combine the two, before ascertaining declension, ascensional difference, etc."[9] Lancelot Wilkinson translated the last of these three verses in a too concise manner to convey the full meaning, and skipped the portion combine the two which the modern Hindi commentary has brought to the fore. According to the Hindi commentary, the final value of period of precession should be obtained by combining +199669 revolutions of ayana with −30000 revolutions of sampaat to get +169669 per Kalpa, i.e. one revolution in 25461 years, which is near the modern value of 25771 years.

Moreover, Munjāla's value gives a period of 21636 years for ayana's motion, which is the modern value of precession when anomalistic precession is also taken into account. The latter has a period of 136000 years now, but Bhāskar-II gives its value at 144000 years (30000 in a Kalpa), calling it sampāt. Bhāskar-II did not give any name of the final term after combining the negative sampāt with the positive ayana. But the value he gave indicates that by ayana he meant precession on account of the combined influence of orbital and anomalistic precessions, and by sampāt he meant the anomalistic period, but defined it as equinox. His language is a bit confused, which he clarified in his own Vāsanābhāshya commentary Siddhānta Shiromani[10] by saying that Suryasiddhanta was not available and he was writing on the basis of hearsay. Bhāskar-II did not give his own opinion, he merely cited Suryasiddhanta, Munjāla and unnamed "others".

Extant Suryasiddhanta supports the notion of trepidation within a range of ±27° at the rate of 54" per year according to traditional commentators, but Burgess opined that the original meaning must have been of a cyclical motion, for which he quoted the Suryasiddhanta mentioned by Bhāskar II.[11]

Yu Xi

Yu Xi (fourth century CE) was the first Chinese astronomer to mention precession. He estimated the rate of precession as 1° in 50 years (Pannekoek 1961, p. 92).

Middle Ages and Renaissance

In medieval Islamic astronomy, the Zij-i Ilkhani compiled at the Maragheh observatory set the precession of the equinoxes at 51 arc seconds per annum, which is very close to the modern value of 50.2 arc seconds.[12]

In the Middle Ages, Islamic and Latin Christian astronomers treated "trepidation" as a motion of the fixed stars to be added to precession. This theory is commonly attributed to the Arab astronomer Thabit ibn Qurra, but the attribution has been contested in modern times. Nicolaus Copernicus published a different account of trepidation in De revolutionibus orbium coelestium (1543). This work makes the first definite reference to precession as the result of a motion of the Earth's axis. Copernicus characterized precession as the third motion of the earth.

Modern period

Over a century later precession was explained in Isaac Newton's Philosophiae Naturalis Principia Mathematica (1687) to be a consequence of gravitation (Evans 1998, p. 246). However, Newton's original precession equations did not work and were revised considerably by Jean le Rond d'Alembert and subsequent scientists.

In 1825, Baron Georges Cuvier, citing work of Jean Baptiste Delambre estimates page 163 the length of a full cycle at 25,960 years, comparing Hipparchus' with measurements taken about 1800 AD. This is less than half of a percent from the current accepted value[13].

Hipparchus' discovery

Hipparchus gave an account of his discovery in On the Displacement of the Solsticial and Equinoctial Points (described in Almagest III.1 and VII.2). He measured the ecliptic longitude of the star Spica during lunar eclipses and found that it was about 6° west of the autumnal equinox. By comparing his own measurements with those of Timocharis of Alexandria (a contemporary of Euclid who worked with Aristillus early in the 3rd century BC), he found that Spica's longitude had decreased by about 2° in about 150 years. He also noticed this motion in other stars. He speculated that only the stars near the zodiac shifted over time. Ptolemy called this his "first hypothesis" (Almagest VII.1), but did not report any later hypothesis Hipparchus might have devised. Hipparchus apparently limited his speculations because he had only a few older observations, which were not very reliable.

Why did Hipparchus need a lunar eclipse to measure the position of a star? The equinoctial points are not marked in the sky, so he needed the Moon as a reference point. Hipparchus had already developed a way to calculate the longitude of the Sun at any moment. A lunar eclipse happens during Full moon, when the Moon is in opposition. At the midpoint of the eclipse, the Moon is precisely 180° from the Sun. Hipparchus is thought to have measured the longitudinal arc separating Spica from the Moon. To this value, he added the calculated longitude of the Sun, plus 180° for the longitude of the Moon. He did the same procedure with Timocharis' data (Evans 1998, p. 251). Observations like these eclipses, incidentally, are the main source of data about when Hipparchus worked, since other biographical information about him is minimal. The lunar eclipses he observed, for instance, took place on April 21, 146 BC, and March 21, 135 BC (Toomer 1984, p. 135 n. 14).

Hipparchus also studied precession in On the Length of the Year. Two kinds of year are relevant to understanding his work. The tropical year is the length of time that the Sun, as viewed from the Earth, takes to return to the same position along the ecliptic (its path among the stars on the celestial sphere). The sidereal year is the length of time that the Sun takes to return to the same position with respect to the stars of the celestial sphere. Precession causes the stars to change their longitude slightly each year, so the sidereal year is longer than the tropical year. Using observations of the equinoxes and solstices, Hipparchus found that the length of the tropical year was 365+1/4−1/300 days, or 365.24667 days (Evans 1998, p. 209). Comparing this with the length of the sidereal year, he calculated that the rate of precession was not less than 1° in a century. From this information, it is possible to calculate that his value for the sidereal year was 365+1/4+1/144 days (Toomer 1978, p. 218). By giving a minimum rate he may have been allowing for errors in observation.

To approximate his tropical year Hipparchus created his own lunisolar calendar by modifying those of Meton and Callippus in On Intercalary Months and Days (now lost), as described by Ptolemy in the Almagest III.1 (Toomer 1984, p. 139). The Babylonian calendar used a cycle of 235 lunar months in 19 years since 499 BC (with only three exceptions before 380 BC), but it did not use a specified number of days. The Metonic cycle (432 BC) assigned 6,940 days to these 19 years producing an average year of 365+1/4+1/76 or 365.26316 days. The Callippic cycle (330 BC) dropped one day from four Metonic cycles (76 years) for an average year of 365+1/4 or 365.25 days. Hipparchus dropped one more day from four Callipic cycles (304 years), creating the Hipparchic cycle with an average year of 365+1/4−1/304 or 365.24671 days, which was close to his tropical year of 365+1/4−1/300 or 365.24667 days. The three Greek cycles were never used to regulate any civil calendar—they only appear in the Almagest in an astronomical context.

We find Hipparchus' mathematical signatures in the Antikythera Mechanism, an ancient astronomical computer of the 2nd century BC. The mechanism is based on a solar year, the Metonic Cycle, which is the period the Moon reappears in the same star in the sky with the same phase (full Moon appears at the same position in the sky approximately in 19 years), the Callipic cycle (which is four Metonic cycles and more accurate), the Saros cycle and the Exeligmos cycles (three Saros cycles for the accurate eclipse prediction). The study of the Antikythera Mechanism proves that the ancients have been using very accurate calendars based on all the aspects of solar and lunar motion in the sky. In fact the Lunar Mechanism which is part of the Antikythera Mechanism depicts the motion of the Moon and its phase, for a given time, using a train of four gears with a pin and slot device which gives a variable lunar velocity that is very close to the second law of Kepler, i.e. it takes into account the fast motion of the Moon at perigee and slower motion at apogee. This discovery proves that Hipparchus mathematics were much more advanced than Ptolemy describes in his books, as it is evident that he developed a good approximation of Kepler΄s second law.

Mithraic question

Mithraism was a mystery religion or school based on the worship of the god Mithras. Many underground temples were built in the Roman Empire from about the 1st century BC to the 5th century AD. Understanding Mithraism has been made difficult by the near-total lack of written descriptions or scripture; the teachings must be reconstructed from iconography found in mithraea (a mithraeum was a cave or underground meeting place that often contained bas reliefs of Mithras, the zodiac and associated symbols). Until the 1970s most scholars followed Franz Cumont in identifying Mithras with the Persian god Mithra. Cumont's thesis was re-examined in 1971, and Mithras is now believed to be a syncretic deity only slightly influenced by Persian religion.

Mithraism is recognized as having pronounced astrological elements, but the details are debated. One scholar of Mithraism, David Ulansey, has interpreted Mithras (Mithras Sol Invictus - the unconquerable sun) as a second sun or star that is responsible for precession. He suggests the cult may have been inspired by Hipparchus' discovery of precession. Part of his analysis is based on the tauroctony, an image of Mithras sacrificing a bull, found in most of the temples. According to Ulansey, the tauroctony is a star chart. Mithras is a second sun or hyper-cosmic sun and or a constellation Perseus, and the bull is Taurus, a constellation of the zodiac. In an earlier astrological age, the vernal equinox had taken place when the Sun was in Taurus. The tauroctony, by this reasoning, commemorated Mithras-Perseus ending the "Age of Taurus" (about 2000 BC based on the Vernal Equinox - or about 11,500 BC based on the Autumnal Equinox).

The iconography also contains two torch bearing boys (Cautes and Cautopates) on each side of the zodiac. Ulansey, and Walter Cruttenden in his book Lost Star of Myth and Time, interpret these to mean ages of growth and decay, or enlightenment and darkness; primal elements of the cosmic progression. Thus Mithraism is thought to have something to do with the changing ages within the precession cycle or Great Year (Plato's term for one complete precession of the equinox).

Changing pole stars

Precession of Earth's axis around the north ecliptical pole
Precession of Earth's axis around the south ecliptical pole

A consequence of the precession is a changing pole star. Currently Polaris is extremely well-suited to mark the position of the north celestial pole, as Polaris is a moderately bright star with a visual magnitude of 2.1 (variable), and it is located within a half degree of the pole.[1]

On the other hand, Thuban in the constellation Draco, which was the pole star in 3000 BC, is much less conspicuous at magnitude 3.67 (one-fifth as bright as Polaris); today it is invisible in light-polluted urban skies.

The brilliant Vega in the constellation Lyra is often touted as the best north star (it fulfilled that role around 12,000 BC and will do so again around the year 14,000); however, it never comes closer than 5° to the pole.

When Polaris becomes the north star again around 27,800, due to its proper motion it then will be farther away from the pole than it is now, while in 23,600 BC it came closer to the pole.

It is more difficult to find the south celestial pole in the sky at this moment, as that area is a particularly bland portion of the sky, and the nominal south pole star is Sigma Octantis, which with magnitude 5.5 is barely visible to the naked eye even under ideal conditions. That will change from the 80th to the 90th centuries, however, when the south celestial pole travels through the False Cross.

This situation also is seen on a star map. The orientation of the south pole is moving toward the Southern Cross constellation. For the last 2,000 years or so, the Southern Cross has nicely pointed to the south pole. By consequence, the constellation is no longer visible from subtropical northern latitudes, as it was in the time of the ancient Greeks.

Polar shift and equinoxes shift

Precessional movement as seen from 'outside' the celestial sphere
Same picture as above, but now from (near) Earth perspective

The figures to the right attempt to explain the relation between the precession of the Earth's axis and the shift in the equinoxes. These figures show the position of the Earth's axis on the celestial sphere, a fictitious sphere which places the stars according to their position as seen from Earth, regardless of their actual distance. The first image shows the celestial sphere from the outside, with the constellations in mirror image. The second figure shows the perspective of a near-Earth position as seen through a very wide angle lens (from which the apparent distortion arises).

The rotation axis of the Earth describes, over a period of 25,700 years, a small circle (blue) among the stars, centered on the ecliptic north pole (the blue E) and with an angular radius of about 23.4°, an angle known as the obliquity of the ecliptic. The direction of precession is opposite to the daily rotation of the Earth on its axis. The orange axis was the Earth's rotation axis 5,000 years ago, when it pointed to the star Thuban. The yellow axis, pointing to Polaris, marks the axis now.

The equinoxes occur where the celestial equator intersects the ecliptic (red line), that is, where the Earth's axis is perpendicular to the line connecting the centers of the Sun and Earth. (Note that the term "equinox" here refers to a point on the celestial sphere so defined, rather than the moment in time when the Sun is overhead at the Equator, though the two meanings are related.) When the axis precesses from one orientation to another, the equatorial plane of the Earth (indicated by the circular grid around the equator) moves. The celestial equator is just the Earth's equator projected onto the celestial sphere, so it moves as the Earth's equatorial plane moves, and the intersection with the ecliptic moves with it. The positions of the poles and equator on Earth do not change, only the orientation of the Earth against the fixed stars.

As seen from the orange grid, 5,000 years ago, the vernal equinox was close to the star Aldebaran of Taurus. Now, as seen from the yellow grid, it has shifted (indicated by the red arrow) to somewhere in the constellation of Pisces.

Still pictures like these are only first approximations as they do not take into account the variable speed of the precession, the variable obliquity of the ecliptic, the planetary precession (which is a slow rotation of the ecliptic plane itself, presently around an axis located on the plane, with longitude 174°.8764) and the proper motions of the stars.

Diagram showing the westward shift of the vernal equinox among the stars over the past six millennia


Cause

The precession of the equinoxes is caused by the gravitational forces of the Sun and the Moon, and to a lesser extent other bodies, on the Earth. It was first explained by Sir Isaac Newton.[14]

Axial precession is similar to the precession of a spinning top. In both cases, the applied force is due to gravity. For a spinning top, this force tends to be almost parallel to the rotation axis. For the Earth, however, the applied forces of the Sun and the Moon are nearly perpendicular to the axis of rotation.

The Earth is not a perfect sphere but an oblate spheroid, with an equatorial diameter about 43 kilometers larger than its polar diameter. Because of the Earth's axial tilt, during most of the year the half of this bulge that is closest to the Sun is off-center, either to the north or to the south, and the far half is off-center on the opposite side. The gravitational pull on the closer half is stronger, since gravity decreases with distance, so this creates a small torque on the Earth as the Sun pulls harder on one side of the Earth than the other. The axis of this torque is roughly perpendicular to the axis of the Earth's rotation so the axis of rotation precesses. If the Earth were a perfect sphere, there would be no precession.

This average torque is perpendicular to the direction in which the rotation axis is tilted away from the ecliptic pole, so that it does not change the axial tilt itself. The magnitude of the torque from the sun (or the moon) varies with the gravitational object's alignment with the Earth's spin axis and approaches zero when it is orthogonal.

Although the above explanation involved the Sun, the same explanation holds true for any object moving around the Earth, along or close to the ecliptic, notably, the Moon. The combined action of the Sun and the Moon is called the lunisolar precession. In addition to the steady progressive motion (resulting in a full circle in about 25,700 years) the Sun and Moon also cause small periodic variations, due to their changing positions. These oscillations, in both precessional speed and axial tilt, are known as the nutation. The most important term has a period of 18.6 years and an amplitude of less than 20 seconds of arc.

In addition to lunisolar precession, the actions of the other planets of the solar system cause the whole ecliptic to rotate slowly around an axis which has an ecliptic longitude of about 174° measured on the instantaneous ecliptic. This so-called planetary precession shift amounts to a rotation of the ecliptic plane of 0.47 seconds of arc per year (more than a hundred times smaller than lunisolar precession). The sum of the two precessions is known as the general precession.

Equations

Tidal force on Earth due to the Sun, Moon, or a planet

The tidal force on Earth due a perturbing body (Sun, Moon or planet) is the result of the inverse-square law of gravity, whereby the gravitational force of the perturbing body on the side of Earth nearest it is greater than the gravitational force on the far side. If the gravitational force of the perturbing body at the center of Earth (equal to the centrifugal force) is subtracted from the gravitational force of the perturbing body everywhere on the surface of Earth, only the tidal force remains. For precession, this tidal force takes the form of two forces which only act on the equatorial bulge outside of a pole-to-pole sphere. This couple can be decomposed into two pairs of components, one pair parallel to Earth's equatorial plane toward and away from the perturbing body which cancel each other, and another pair parallel to Earth's rotational axis, both toward the ecliptic plane.[15] The latter pair of forces creates the following torque vector on Earth's equatorial bulge:[3]

\overrightarrow{T} = \frac{3Gm}{r^3}(C-A)\sin\delta\cos\delta\begin{pmatrix}\sin\alpha\\-\cos\alpha\\0\end{pmatrix}

where

Gm = standard gravitational parameter of the perturbing body
r = geocentric distance to the perturbing body
C = moment of inertia around Earth's axis of rotation
A = moment of inertia around any equatorial diameter of Earth
C−A = moment of inertia of Earth's equatorial bulge (C>A)
δ = declination of the perturbing body (north or south of equator)
α = right ascension of the perturbing body (east from vernal equinox)

The three unit vectors of the torque at the center of the Earth (top to bottom) are x on a line within the ecliptic plane (the intersection of Earth's equatorial plane with the ecliptic plane) directed toward the vernal equinox, y on a line in the ecliptic plane directed toward the summer solstice (90° east of x), and z on a line directed toward the north pole of the ecliptic.

The value of the three sinusoidal terms in the direction of x (sinδ cosδ sinα) for the Sun is a sine squared waveform varying from zero at the equinoxes (0°, 180°) to 0.36495 at the solstices (90°, 270°). The value in the direction of y (sinδ cosδ (−cosα)) for the Sun is a sine wave varying from zero at the four equinoxes and solstices to ±0.19364 (slightly more than half of the sine squared peak) halfway between each equinox and solstice with peaks slightly skewed toward the equinoxes (43.37°(−), 136.63°(+), 223.37°(−), 316.63°(+)). Both solar waveforms have about the same peak-to-peak amplitude and the same period, half of a revolution or half of a year. The value in the direction of z is zero.

The average torque of the sine wave in the direction of y is zero for the Sun or Moon, so this component of the torque does not affect precession. The average torque of the sine squared waveform in the direction of x for the Sun or Moon is:

T_x = \frac{3}{2}\frac{Gm}{a^3 (1-e^2)^{3/2}}(C-A)\sin\epsilon\cos\epsilon

where

a = semimajor axis of Earth's (Sun's) orbit or Moon's orbit
e = eccentricity of Earth's (Sun's) orbit or Moon's orbit

and 1/2 accounts for the average of the sine squared waveform, a^3 (1-e^2)^{3/2} accounts for the average distance cubed of the Sun or Moon from Earth over the entire elliptical orbit,[16] and \epsilon\,\! (the angle between the equatorial plane and the ecliptic plane) is the maximum value of δ for the Sun and the average maximum value for the Moon over an entire 18.6 year cycle.

Precession is:

\frac{d\psi}{dt} = \frac{T_x}{C\omega\sin\epsilon}

where ω is Earth's angular velocity and Cω is Earth's angular momentum. Thus the first order component of precession due to the Sun is:[3]

\frac{d\psi_S}{dt} = \frac{3}{2}\left[\frac{Gm}{a^3 (1-e^2)^{3/2}}\right]_S\left[\frac{(C-A)}{C}\frac{\cos\epsilon}{\omega}\right]_E

whereas that due to the Moon is:

\frac{d\psi_L}{dt} = \frac{3}{2}\left[\frac{Gm(1-1.5\sin^2 i)}{a^3 (1-e^2)^{3/2}}\right]_L\left[\frac{(C-A)}{C}\frac{\cos\epsilon}{\omega}\right]_E

where i is the angle between the plane of the Moon's orbit and the ecliptic plane. In these two equations, the Sun's parameters are within square brackets labeled S, the Moon's parameters are within square brackets labeled L, and the Earth's parameters are within square brackets labeled E. The term (1-1.5\sin^2 i) accounts for the inclination of the Moon's orbit relative to the ecliptic. The term (C−A)/C is Earth's dynamical ellipticity or flattening, which is adjusted to the observed precession because Earth's internal structure is not known with sufficient detail. If Earth were homogeneous the term would equal its third eccentricity squared,[17]

e''^2 = \frac{\mathrm{a}^2 - \mathrm{c}^2}{\mathrm{a}^2 + \mathrm{c}^2}

where a is the equatorial radius (6378137 m) and c is the polar radius (6356752 m), so e'' ²=0.003358481.

Applicable parameters for J2000.0 rounded to seven significant digits (excluding leading 1) are:[18][19]

Sun Moon Earth
Gm=1.3271244 × 1020 m³/s² Gm=4.902799 × 1012 m³/s² (C−A)/C=0.003273763
a=1.4959802 × 1011 m a=3.833978 × 108 m ω=7.292115 × 10−5 rad/s
e=0.016708634 e=0.05554553 \epsilon\,\!=23.43928°
i= 5.156690°

which yield

dψS/dt = 2.450183 × 10−12 /s
dψL/dt = 5.334529 × 10−12 /s

both of which must be converted to "/a (arcseconds/annum) by the number of arcseconds in 2π radians (1.296 × 106"/2π) and the number of seconds in one annum (a Julian year) (3.15576 × 107s/a):

dψS/dt = 15.948788"/a   vs   15.948870"/a from Williams[3]
dψL/dt = 34.723638"/a   vs   34.457698"/a from Williams

The solar equation is a good representation of precession due the Sun because Earth's orbit is close to an ellipse, being only slightly perturbed by the other planets. The lunar equation is not as good a representation of precession due to the Moon because its orbit is greatly distorted by the Sun.

Values

Simon Newcomb's calculation at the end of the nineteenth century for general precession (known as p) in longitude gave a value of 5,025.64 arcseconds per tropical century, and was the generally accepted value until artificial satellites delivered more accurate observations and electronic computers allowed more elaborate models to be calculated. Lieske developed an updated theory in 1976, where p equals 5,029.0966 arcseconds per Julian century. Modern techniques such as VLBI and LLR allowed further refinements, and the International Astronomical Union adopted a new constant value in 2000, and new computation methods and polynomial expressions in 2003 and 2006; the accumulated precession is:

pA = 5,028.796195×T + 1.1054348×T² + higher order terms,

in arcseconds per Julian century, with T, the time in Julian centuries (that is, 36,525 days) since the epoch of 2000.

The rate of precession is the derivative of that:

p = 5,028.796195 + 2.2108696×T + higher order terms

The constant term of this speed corresponds to one full precession circle in 25,772 years.

The precession rate is not a constant, but is (at the moment) slowly increasing over time, as indicated by the linear (and higher order) terms in T. In any case it must be stressed that this formula is only valid over a limited time period. It is clear that if T gets large enough (far in the future or far in the past), the T² term will dominate and p will go to very large values. In reality, more elaborate calculations on the numerical model of solar system show that the precessional constants have a period of about 41,000 years, the same as the obliquity of the ecliptic. Note that the constants mentioned here are the linear and all higher terms of the formula above, not the precession itself. That is,

p = A + BT + CT² + …

is an approximation of

p = A + Bsin (2πT/P), where P is the 410-century period.

Theoretical models may calculate the proper constants (coefficients) corresponding to the higher powers of T, but since it is impossible for a (finite) polynomial to match a periodic function over all numbers, the error in all such approximations will grow without bound as T increases. In that respect, the International Astronomical Union chose the best-developed available theory. For up to a few centuries in the past and the future, all formulas do not diverge very much. For up to a few thousand years in the past and the future, most agree to some accuracy. For eras farther out, discrepancies become too large — the exact rate and period of precession may not be computed, even for a single whole precession period.

The precession of Earth's axis is a very slow effect, but at the level of accuracy at which astronomers work, it does need to be taken into account on a daily basis. Note that although the precession and the tilt of Earth's axis (the obliquity of the ecliptic) are calculated from the same theory and thus, are related to each other, the two movements act independently of each other, moving in mutually perpendicular directions.

Precession exhibits a secular decrease due to tidal dissipation from 59"/a to 45"/a (a = annum = Julian year) during the 500 million year period centered on the present. After short-term fluctuations (tens of thousands of years) are averaged out, the long-term trend can be approximated by the following polynomials for negative and positive time from the present in "/a, where T is in billions of Julian years (Ga):[20]

p− = 50.475838 − 26.368583T + 21.890862T²
p+ = 50.475838 − 27.000654T + 15.603265T²

Precession will be greater than p+ by the small amount of +0.135052"/a between +30 Ma and +130 Ma. The jump to this excess over p+ will occur in only 20 Ma beginning now because the secular decrease in precession is beginning to cross a resonance in Earth's orbit caused by the other planets.

According to Ward, when, in about 1,500 million years, the distance of the Moon, which is continuously increasing from tidal effects, has increased from the current 60.3 to approximately 66.5 Earth radii, resonances from planetary effects will push precession to 49,000 years at first, and then, when the Moon reaches 68 Earth radii in about 2,000 million years, to 69,000 years. This will be associated with wild swings in the obliquity of the ecliptic as well. Ward, however, used the abnormally large modern value for tidal dissipation. Using the 620-million year average provided by tidal rhythmites of about half the modern value, these resonances will not be reached until about 3,000 and 4,000 million years, respectively. Long before that time (about 2,100 million years from now), due to the gradually increasing luminosity of the Sun, the oceans of the Earth will have vaporized, which will reduce tidal effects significantly.

See also

References

  1. ↑ 1.0 1.1 1.2 1.3 Astro 101 - Precession of the Equinox, Western Washington University Planetarium, accessed 30 December 2008
  2. ↑ Robert Main, Practical and Spherical Astronomy (Cambridge: 1863) pp.203–4.
  3. ↑ 3.0 3.1 3.2 3.3 James G. Williams, "Contributions to the Earth's obliquity rate, precession, and nutation", Astronomical Journal 108 (1994) 711–724, pp.712&716. All equations are from Williams.
  4. ↑ IAU 2006 Resolution B1: Adoption of the P03 Precession Theory and Definition of the Ecliptic
  5. ↑ Dennis Rawlins, Continued fraction decipherment: the Aristarchan ancestry of Hipparchos' yearlength & precession DIO (1999) 30–42.
  6. ↑ Neugebauer, O. "The Alleged Babylonian Discovery of the Precession of the Equinoxes", Journal of the American Oriental Society, Vol. 70, No. 1. (Jan. – Mar., 1950), pp. 1–8.
  7. ↑ Milbrath, S. "Just How Precise is Maya Astronomy?", Institute of Maya Studies newsletter, December 2007.
  8. ↑ Siddhānta-shiromani, Golādhyāya, section-VI, verses 17-19
  9. ↑ Translation of the Surya Siddhānta by Pundit Bāpu Deva Sāstri and of the Siddhānta Siromani by the Late Lancelot Wilkinson revised by Pundit Bāpu Deva Sāstri, printed by C B Lewis at Baptist Mission Press, Calcutta, 1861; Siddhānta Shiromani Hindi commentary by Pt Satyadeva Sharmā, Chowkhambā Surbhārati Prakāshan, Varanasi, India.
  10. ↑ Vāsanābhāshya commentary Siddhānta Shiromani (published by Chowkhamba)
  11. ↑ cf. Suryasiddhanta, commentary by E. Burgess, ch.iii, verses 9-12.
  12. ↑ Rufus, W. C. (May 1939), "The Influence of Islamic Astronomy in Europe and the Far East", Popular Astronomy 47 (5): 233–238 [236] 
  13. ↑ Ancient records of astronomy are discussed in Baron Georges Cuvier A Discourse on the Revolutions of the Surface of the Globe (1825) in the chapter titled "The astronomical monuments of the ancients." Part of the discussion in the chapter concerns the date of the Dendera zodiac, which Cuvier suggested is 123 AD to 147 AD page 170 and page 172. The current consensus is around 50 BC.
  14. ↑ The Columbia Electronic Encyclopedia, 6th ed., 2007
  15. ↑ Ivan I. Mueller, Spherical and practical astronomy as applied to geodesy (New York: Frederick Unger, 1969) 59.
  16. ↑ G. Boué & J. Laskar, "Precession of a planet with a satellite", Icarus 185 (2006) 312–330, p.329.
  17. ↑ George Biddel Airy, Mathematical tracts on the lunar and planetary theories, the figure of the earth, precession and nutation, the calculus of variations, and the undulatory theory of optics (third edititon, 1842) 200.
  18. ↑ J.L. Simon et al., "Numerical expressions for precession formulae and mean elements for the Moon and the planets", Astronomy and Astrophyics 282 (1994) 663-683.
  19. ↑ Dennis D. McCarthy, IERS Technical Note 13 – IERS Standards (1992) (Postscript, use PS2PDF).
  20. ↑ J. Laskar et al., "A long-term numerical solution for the insolation quantities of the Earth", Astronomy and Astrophysics 428 (2004) 261-285, pp.276 & 278.

External links